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The gauge-boson sector of perturbative scalar electrodynamics is investi- 
gated in detail as a testing ground for resummation methods in hot gauge theo- 
ries. It also serves as a simple non-trivial reference system for the non-Abelian 
gluon plasma. The complete next-to-leading order contributions to the polar- 
ization tensor are obtained within the resummation scheme of Braaten and 
Pisarski. The simpler scheme proposed recently by Arnold and Espinosa is 
shown to apply to static quantities only, whereas Braaten-Pisarski resum- 
mation turns out to need modification for collective phenomena close to the 
light-cone. Finally, a recently proposed resummation of quasi-particle damp- 
ing contributions is assessed critically. 



1. Introduction 



Perturbative thermal gauge field theories at ultra-relativistic temperatures |T| are 
hoped to be applicable for the description of phenomena associated with the envisaged 
formation of a quark-gluon plasma in high-energy ion collisions and for the physics of the 
early universe. In recent years, a lot of work has been put into developing a systematic 
perturbation scheme, which has turned out to require a resummation of the usual series 
of loop diagrams. 

In non-Abelian gauge theories it was found by Linde that there exists a barrier 
raised by infrared divergences associated with unscreened static magnetic fields, but even 
before this barrier is hit, conventional perturbation theory breaks down for certain mass 
scales much smaller than the temperature. In hot quantum chromodynamics (QCD), 
which is the fundamental theory relevant for the hypothetical quark-gluon plasma, it 
is for example the spectrum of quasi-particles whose perturbative treatment requires a 
resummed perturbation theory. After some failed attempts to calculate dissipative proper- 
ties of gluonic plasma excitations , such a resummation scheme was developed finally by 
Braaten and Pisarski Q. Also in spontaneously broken theories, there has been a flurry of 
activity centered around the issue of determining corrections to the effective potentials of 
scalar fields which describe the cosmological phase transition from the high-temperature 
symmetric to the low-temperature broken phase. These developments, which have been 
pioneered in the early seventies ||, are summarized e.g. in Ref. ||. 

The need for a resummation of the ordinary loop expansion is not peculiar to non- 
Abelian theories, but can be studied also in the Abelian case. One of the simplest models 
is scalar electrodynamics. As Higgs model, it has already been used extensively to inves- 
tigate the nature of the phase transition at finite temperature, e.g. in Ref. f7|. There the 
scalar sector is of central interest. In this paper we shall concern ourselves with the sector 
of the gauge fields, and we shall use hot scalar electrodynamics to study en miniature 
the issues that have been raised in the case of hot QCD. Scalar electrodynamics appears 
to be a potentially interesting toy model for purely gluonic QCD, since it involves self- 
interacting bosons and a gauge-field sector. In hot QCD, the resummed perturbation 
theory of Ref. M has been employed to derive a consistent, gauge-independent damping 
constant for the lowest-energy plasmon mode || as well as for energetic quasi-particles 
0, and also some next-to-leading order results: corrections to the plasma frequency |L0| 
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and to Debye screening ]TT[ . In the much simpler case of hot scalar electrodynamics, we 
shall be able to give a fairly complete picture for the photonic quasi-particle spectrum. 



In doing so, we shall largely follow the resummation program of Braaten and Pisarski 
0]. After presenting the leading-order results in Sect. 2, this scheme is reviewed for the 
simple case of scalar electrodynamics in Sect. 3. In Sect. 4 we first specialize to the 
static limit and obtain next-to-leading order results for Debye screening and magnetic 
permeability. Here the simpler resummation scheme of Ref. |J, which has originally 
been put forward for the resummation of the static effective potential, agrees with full 
resummation. In Sect. 5, the limit of long- wavelength oscillations is considered and the 
next-to-leading order result for the plasma frequency is found. In this case, a naive 
application of the resummation scheme of Arnold and Espinosa is shown to fail, because 
nonstatic modes can no longer be neglected. We also show to what extent classical 
considerations can explain the value of the plasma frequency obtained before. In Sect. 6 
we exhibit the complete next-to-leading order results for the polarization tensor, deriving 
the corrected spectrum of propagating photonic quasi-particles. Here we find a qualitative 
change in the case of the longitudinal plasmons, which turn out to exist only for a finite 
range of frequencies or momenta. The canonical resummation according to Braaten and 
Pisarski actually breaks down in this example due to singularities at the light-cone, but 
can be amended so that this upper bound for the frequencies of longitudinal plasmons can 
be calculated accurately. A similar result is shown to hold in the case of QCD. Finally, in 
Sect. 7 we consider the additional resummation of scalar damping contributions, critically 
assessing a recent proposal in Ref. O]. Sect. 8 contains our conclusions. 



2. Leading-order results 



Scalar electrodynamics with a scalar potential leading to spontaneous symmetry break- 
ing has been extensively used as a toy model to study symmetry restoration at high tem- 
perature and the corresponding phase transition which occurs when the temperature is 
lowered 0. In this paper we shall focus our attention to the gauge field sector of this 
model rather than the effective potential of the scalars. For this we shall for simplicity 
consider charged scalar particles without self-interactions. At sufficiently high tempera- 
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ture, the bare mass of the scalar particles can be neglected, so our starting point is the 



Lagrangian [13 



C = (D^)* D»<p - - F, U F^ - — {d»A»f , (2.1) 

with the covariant derivative = + ieA^ and the Abelian field tensor = dy,A v — 

duAy,. The signature of the Minkowski metric g^ v is H . In fl2.1|) we have included a 

gauge breaking term corresponding to general covariant gauges; the ghost term is omitted 
since it decouples from the rest. 

We shall mainly use the Matsubara formalism, where Green's functions are first de- 
fined for discrete imaginary frequencies proportional to 2niT. After the evaluation of all 
frequency sums, they are finally extended to real external frequencies by an appropri- 
ate analytic continuation. In the first chapters we shall have to consider exclusively the 
continuation of two-point functions, for which we choose retarded boundary conditions 
prescribing Q = 2ninT — > uj + is [0, @ • 



The theory (|2.1| ) has two propagators and two vertices. Writing a four-momentum 
as Q = (iui n ,q), u) n = 27inT, the bare scalar propagator is 5° = —1/Q 2 , and the bare 
photon propagator reads G® = g^u/Q 2 + (« — 1)D^ V /Q 2 . Here D is the fourth matrix of 
the following basis for symmetric Lorentz tensors, 

VoV QoV+VoQ QoQ , 

A = „-B-D , B = — , ^ , D^yi , (2.2) 

where V = Q 2 U — (UQ) Q is a projector to longitudinal fields built from the rest-frame 
velocity of the heat bath U, which in the following we choose as U — (1,0). A and B are 
transverse with respect to Q M ; A is transverse also with respect to the three- momentum 
q. The translation into the notation of the textbook of Kapusta U is given by A = —Pt 
and B = —Pl- The vertex which couples a scalar line (Q ingoing, P outgoing) with one 
photon is — e(Q + P) M , and that connecting to two photon lines is 2e 2 g fiu with the Lorentz 
indices carried by the two photons. 

At one-loop level, the propagators are given by Dyson's equation 

G = G° + G°nG and S = S° + S° E S , (2.3) 

where 77 and S are the self energies of photon and scalar, respectively. The corresponding 
Feynman diagrams are depicted in Fig. 1. 
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Figure 1: Propagators dressed by hard thermal loops. Dotted lines represent photons, and solid lines 
stand for scalar particles. 

The bare 1-loop expression for the photon self-energy reads 

n^(Q) = 6 -mV» + e 2 E K2 {K _ Q) 2 with m =Y > (2-4) 

where the blank sum symbol stands for (2Tr)~ 3 Jd 3 k T^ n - We shall always use K as 
loop-momentum variable. In contrast to the self-energy of gluons in QCD, ( |2.4j ) is gauge 
parameter independent, and also somewhat simpler. i7 Atl/ is transverse, Q fl II fll/ = which 
in Abelian theories holds for any linear gauge choice (in non-Abelian theories the finite- 
temperature gluon self-energy is, in general, non-transverse even in covariant gauges). It 
therefore decomposes into 

n v» = A ^jj t + B ^jj £ (2.5) 
with n t = Tr An/2 and i7 £ = Tr BII. 

In the high-temperature limit, the leading order contributions are given by 

„ t = -%n m = w (i - f) (i - 9l ln «i±l) , (2 . 6) 

T V T J \ 2q Q - qj 

n t = i(3m 2 -il/) . (2.7) 

For the comparison with subleading results to be discussed later on, in Fig. 2 the 
real part of n e is displayed as a function of real Q = uj and q (according to Q2.7| ) II t is 
simply given by the inverted picture, appropriately relabelled). Notice in particular the 
non-analytic behavior at the origin, where the value of lie (and also of n t ) depends on 
whether it is approached from space-like, time-like, or light-like directions. 

The scalar self-energy on the other hand is given by 

- 2 2^2^/ 3 -a 2(a-l)KQ\ eT 

a = " S [ k^K-QY + K 4 (K — Q) 2 ) mth " s T ■ (2 ' 8) 
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Figure 2: The real part of Z7^(cl>, q) in units of e 2 T 2 . Also given is the intersection with the surface 

2 2 

— which gives the dispersion curve for longitudinal plasmons. 

(Including <y2 4 -self-interactions by £ — > £ — j(ip*ip) 2 would only change fi 2 — > e 2 T 2 /4 + 
AT 2 /12.) In contrast to the photon self-energy, S does depend on the gauge parameter 
a, but the only contribution oc T 2 is the first term, which is again gauge- independent. 

Including only the leading order contributions oc e 2 T 2 , the scalar propagator according 
to ( |2.3| ) is simply 

with a constant (thermal) mass term. The photon propagator, however, has a richer 
structure, and reads in the tensorial basis (12.21) 



G ^> = a ^q^WT) + ^Q^nW) + D ^ • (210) 

with II t and II \ at leading order as given by (|2.6| , |2~7| ) . 

The ultra-relativistic photon propagator (p.lCQ has been first derived within classical 
kinetic theory by Silin and within field theory by Fradkin |16[|. Exactly the same 



expressions were later found in high-temperature QCD |T7| , 18 1, the only difference be- 
ing a replacement e 2 — > g 2 (N + Nf/2) for SU(iV) with Nf flavors. There are now two 
physical structure functions in the gauge-field propagator. The one associated with the 
spatially transverse tensor A corresponds to transverse photons (or gluons), whereas the 
one associated with B describes a new collective mode, the so-called plasmon mode. The 
poles of these structure functions correspond to the normal modes of the gauge field sector 
of the ultra-relativistic plasma. For real frequencies and momenta, these give the spec- 
trum of propagating photonic quasi-particles. It starts at the minimal (plasma) frequency 
uj(q = 0) = eT/3 = m. For larger frequencies the two branches differ. At q — > oo the spa- 



tially transverse branch approaches the form yq 2 + rn 2 ^ with an asymptotic thermal mass 
= eT/v^ 1 , whereas the longitudinal branch approaches the light-cone exponentially. 
This longitudinal branch is also shown in Fig. 2. It disappears from the spectrum even- 
tually because with increasing momentum the residue of the corresponding pole in ( 2.10Q 



vanishes exponentially. Perturbing the plasma with frequencies u < m does not give rise 
to propagating quasi-particles, but results in dynamical screening, both in the spatially 
transverse and in the longitudinal mode. In the static limit, there is screening only in the 
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longitudinal sector corresponding to Debye screening of longitudinal electric fields with 
screening mass m e \, = eT/y/S 1 , whereas static magnetic fields remain unscreened. 

The main theme of this work will be to determine this spectrum of the photonic 
quasi-particles and the screening properties beyond the level of the leading temperature 
contributions. 



3. Resummation of hard thermal loops 



In perturbative quantum field theories it sometimes happens that higher-order loop 
diagrams are not suppressed by correspondingly high powers of the coupling constants, 
and that therefore the perturbation series has to be reorganized in order to actually be 
perturbative. In thermal field theories this is a ubiquitious phenomenon. In particular 
the zero modes of massless Bose particles are sensitive to the appearance of thermal 
masses. Repeated self-energy insertions cause higher-order diagrams to become more and 
more infra-red singular, but summing them up according to Fig. 1 the resulting dressed 
propagators are well-behaved. In the case of the thermodynamic potential, this particular 
resummation is known as "ring resummation" and dates back to the work of Gell-Mann 



and Brueckner 19 



In the case of the spectrum of quasi-particle excitations, the necessity of a resummed 
perturbation theory has become apparent in gauge theories by the failure of the unim- 
proved perturbation theory to give gauge-independent results for the damping constants 
of thermal quasi-particle excitations . A systematic perturbation theory was developed 
most notably by Braaten and Pisarski Q. These authors have shown that Green's func- 
tions involving soft external momenta (i.e. of the same order of magnitude as the plasma 
frequency) require the resummation of all so-called hard thermal loops (HTL's) |20| [|. 
The latter are the contributions from one-loop diagrams that are dominated by hard loop 
momenta ~ T. In renormalizable theories they go like T 2 for T — > oo. They can compen- 
sate for powers of the coupling as follows. Consider in particular the dressed propagators 
as given by Fig. 1. In the diagrams on the right-hand-side of Fig. 1 the self-energy part 
contributes a factor ~ e 2 T 2 . If the bare propagator attached to it carries momentum 
~ eT, it brings in a factor e~ 2 T~ 2 and the contribution involving the self-energy inser- 
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tion is of relative order 1. Only for larger external momentum the perturbative nature 
of Fig. 1 is restored, and loop-corrections are suppressed by factors of e. Consequently, 
contributions from loop integrals (in any higher order) that probe these soft scales are 
not perturbatively stable, unless (at least) these hard thermal loops are resummed. The 
hard thermal loops themselves are stable, because they are sensitive to hard scales only. 

In non-Abelian gauge theories, the same situation occurs for vertex functions, because 
there the one- loop vertices receive also contributions ~ T 2 , such that HTL vertices are 
of the same order of magnitude as bare ones for soft external momenta. Therefore, they 
have to be included in the resummation as well. 

In scalar electrodynamics, there are no HTL vertices. Superficial power counting M 
rules out (as in QCD) all candidates for HTL vertices, which contain a bare 4-vertex or 
different kinds of lines in the loop. Moreover, scalar loops with an odd number of external 
lines are easily shown to vanish, because they change sign by reversing momentum flows. 
For the remaining scalar loops with an even number n of external photons with momenta 
Qj one arrives at 

^ ^^(^ + perm, of Q,,...^ . (3.1) 

with Q\ + Q 2 + . . . + Q n = 0. The power-counting arguments of Ref. || would estimate 
that ( |3.1D goes like T 2 for T — > oo, so that with Qj ~ eT ( |5TT] ) would become of the order 
e n ~ 2 \Q\ 4: ~ n , the same order of magnitude as bare vertices are on dimensional grounds. 
However, for n — 4 it has been shown in Ref. |^0| that the leading-order contributions 



in ( ^.1| ) in fact cancel, and in Ref. [2T| it was proved by induction that this holds for all 
n > 4 as well. 

In non-Abelian gauge theories with or without fermions, the HTL vertices are essential 
to guarantee the gauge invariance of the generating functional of all HTL's. An explicit 



effective action S e ff has been constructed by Taylor and Wong |22] by starting from its 

(2) 

bilinear part S^ s as determined by the HTL self-energy diagrams and determining higher 
orders in the fields by gauge invariance. A nice representation of this can be found in 



(2) 

Ref. [KJ. In scalar electrodynamics it turns out that S^ s is already gauge invariant so 



(2) 

that S e g = Sg ff , and the uniqueness of this construction, which is presented in some detail 
in App. A, explains the absence of HTL vertices in the Abelian case. 

For soft momenta, the HTL's are of the same order as the tree-level terms, so together 
they represent the 'zeroth order' of the high-temperature limit of the theory. Thus, the 
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effective action, which combines them into a single formula, represents the first term of 
the high-temperatur asymptotics of this theory. In its manifest gauge invariant form, it 
reads 



S cS = S + 6S = f C cS , £ eff = C - mW + \™ 2 jf Y p F Pil — L_ F» X Y> 



(3.2) 



where J 13 = j^dr Jd 3 r and = (1, e) with e 2 = 1. The angular integral J n over the 
directions of e is normalized to one: J n l = 1, J n Y = U . 



When adopting ( |3.2|) as the adequate starting point for a perturbative treatment, one 
must however take care of not to change the underlying theory itself. One is therefore led 
to rewrite S as S = S c ^ — SS, where —5S generates counter-terms that subtract at higher 
loop orders what has been added in at lower ones. 

The first contributions beyond the new zeroth order are contained in 1-loop diagrams, 
but now with resummed internal lines. For the photon self-energy this resummation 
amounts to insert the massive scalar propagators fl2.9|) into the diagrams of Fig.l. The 
general form of the photon Green's function is still given by (|2.10|) by virtue of gauge 
invariance, i.e. transversality. 

The resummed 1PI expressions Hi = TiBII and II t = Tr All/ 2 derive from 

1 . 1 



= -2eV ly E A + 4e2 E^^ A_A > A = 



{K-QY-i* 



and are given more explicitly by 

n e = -2e 2 A + 4e2 E A ~ A 



kg) 



K 



(KQf 

Q 2 



Tig 



where Yl t is determined to the slightly simpler auxiliary quantity 

n g = Ti 9 n = -4e 2 E A + e2 ( V - Q 2 ) E A ~ A 



(3.3) 

(3.4) 
(3.5) 

(3.6) 



In ( |3.3D a term involving Q^Q V has been dropped, because it does not contribute at the 
subleading order under study. 



An alternative version of Hi is 

n e = e 2 (4/i 2 - Q 2 ) J2 A_A + 4e2 E A_ A 



k q] 



(3.7) 
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Special elements of 11^ may be read off from = A^n t + B^ u IIi. In particular, 

n 00 = -^n e , n« = -m t - ^n e = n m - n g . (3.8) 



4. Static screening at next-to-leading order 



In the static limit the spatially-transverse and spatially-longitudinal structure func- 
tions in the photon propagator, as given by (|3~8| ), simplify to 

/7t(0,g) = ~i7«(0,g) , 77,(0, q) = 77 00 (0, q) . (4.1) 

Before resummation, the leading and subleading terms in the high-temperature ex- 
pansion of the static one-loop self-energy (ET4l) read 



iloo(0, q) = + ^- 2 q 2 In °- + . . . = 3m 2 (l + 0(e 2 )) (4.2) 



24tt 2 ^ T 
and 

nu(0, q) = - l -e 2 qT - ^q 2 In °- + . . . = - jjmge (1 + 0(e)) , (4.3) 

where here and in what follows we count orders in e always for soft momenta ~ m cx eT. 
Notice further that the T = part has been assumed to be renormalized at the scale 
a, which contributes a term proportional to \na 2 /Q 2 . This has been combined with a 
similar term proportional to \nQ 2 /T 2 from the temperature-dependent parts. 

The leading contribution in 77oo is the familiar electric (Debye) screening mass m 2 x = 
3m 2 , while there is no screening mass of this order of magnitude in 77^, i.e. no magnetic 
mass. At relative order e there is no contribution in 77 o, but one in 77^ which when 
taken seriously spells trouble. Including the latter in the static transverse propagator 
would give 

which has a space-like pole at q — e 2 T/16. A similar behaviour is found also in QCD, 
sometimes called the Landau ghost of thermal QCD. 
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However, as we have argued in the previous sections, the subleading terms in (|4.2|) 
and ( }4.3| ) are not accurate for q<eT and require resummation of the hard thermal loops. 
Only the terms proportional to e 2 \n(a/T) are to be trusted as they are determined by 
the ultraviolet part of the loop integrals. 



4. A RESUMMATION 



-^oo(0, q) upon resummation of the hard thermal loop, which is just the thermal scalar 
mass [A, reads 

A)o(0,g)=e 2 ]T{4K 2 AA--2A} , (4.5) 

see ( p.4| ), with A and A~ as defined in ( |3.3| ). Evaluating the sum over the Matsubara 
frequencies by means of a contour integral leads to 



#oo(0, q) = — 2 \ dkk : 

7l Z JO 



n[\Jk 2 + /i 2 ' 



kq 



2k + q 



2k -q 



(4.6) 



After separating off the T 2 -contribution, we can write 



#oo(0, q) 



e 2 T 2 e 2 T 



IT 2 JO 



dk 



k 2 + /r 



+ l--ln 

Q 



2k + q 



2k -q 



+ 0(e 2 q 2 \n(T)) . (4.7) 



This can be evaluated by integration by parts, which reveals that the (/-dependence of the 
contribution at order T is completely spurious: 



i7 00 (0, q) = - t^L + 0(e 2 q 2 ln(T)) = 3m 2 (l - A e + (e 2 ] 
3 2n V 47T 



(4- 



The static I7 00 at next-to- leading order is a negative constant, resulting in a decrease 
of the classical value of the electric screening mass. In QCD the corresponding calculation 
has been performed recently by one of the present authors in Ref . [T^j . There the next- 
to- leading order correction to i7oo(0, q) turns out to be a nontrivial function of q, which is 
such that it diverges logarithmically for q 2 — > —m^ , where i7 o(0, q) defines the correction 
term to the Debye mass. Assuming that this divergence is cut-off by a magnetic mass, 
which is expected to arise non-perturbatively at order g 2 T (g being the QCD coupling 
constant), leads to 5m 2 j = 0(g ln(l/g)m 2 ) such that 5m c \ > for small g. Evidently, 
in this quantity there is a qualitative difference between an Abelian and a non-Abelian 
gauge field plasma beyond leading order. 
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In a similar manner, using ( |3.5| ) and ( |3.6|) , one can evaluate the resummed expression 
for n u (0,q), 

n ii (0,q) = e 2 Y,{^k 2 -q 2 )AA- + 6A} . (4.9) 



Its high-temperature limit is found to be 

+ 0(eV) , (4-10) 



which is to be compared with ( |4.3|) . The latter is only accurate for q 3> eT, and in this 
regime it indeed coincides with the resummed result (|4.10| ), since arctan(g/(2|u)) — > n/2 
for qj \i — ► oo. 

For q < eT, the resummed result deviates considerably from the bare one. In particular 
for q — > 0, the former approaches zero like q 2 rather than q, so that there is no longer any 
unphysical pole at space-like momentum. The vanishing of i7jj(0, q —>■ 0) also implies that 
there is no magnetic mass squared of the order of e 3 T 2 . In fact, in the present Abelian 



case it can be shown rigorously that there is no magnetic screening mass [16, 17]. In QCD 



17j the situation is again quite different. Resummation changes the unphysical pole at 
space-like momenta, but does not remove it. However, since this pole arises at momentum 
scale g 2 T, this just points to the relevance of a magnetic mass in the non-Abelian case. 

For finite q ~ eT, 77^(0, q) is always negative, which implies that the magnetic per- 
meability 

(H. lM1 -^ M (4 - n) 

is positive, decreasing monotoneously from e/(127r) at q = to zero for q — > oo. Hence, 
the hot scalar plasma is weakly diamagnetic at distances > l/(eT). 



4. B STATIC RESUMMATION 

Up to now we have strictly followed the resummation scheme outlined by Braaten 
and Pisarski for hot gauge theories (a detailed exposition in the simpler scalar case can be 
found in Ref. [p4|). However, the above calculations can be somewhat simplified by the 
following observations due to Arnold and Espinosa || made in the context of a resummed 
perturbation theory for the finite temperature effective potential in gauge theories. In 
the Matsubara formalism, K 2 has Euclidean signature, — K 2 = (2irnT) 2 + k 2 , so that the 
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momentum can be soft only with n = 0. Accordingly it is sufficient to dress only these 
static modes: 

= wrmj) • (412) 

whereas the correction terms to the nonstatic propagators in 

-f 



(2vmT) 2 + k 
are down by powers of g 2m when 77 ~ g 2 T 2 . 



« (-l) m n m (2ninT,k) 
1 + [(27mT) 2 + P]™ 



(4.13) 



A systematic resummation scheme based on this splitting has been put forward in 
Ref. 0, where the decisive simplifications for gauge theories are due to the fact that the 
hard thermal loop 11(0, k) is a constant mass term. In its lowest order version, it just 
coincides with the well-known ring resummation introduced by Gell-Mann and Brueckner 



0. 



The limitations of this restricted resummation scheme become apparent when one 
considers nonstatic Green's functions. Because of energy conservation at vertices, external 
frequencies are also fixed to Matsubara frequencies, and because of the special treatment 
of the static sector, it is quite impossible to perform an analytic continuation to nonzero 
soft external frequencies in the end. However, this scheme should be sufficient for the 
evaluation of corrections to static quantities like the effective potential pi or screening 
PI , and to Green's functions with hard external frequencies, e.g. damping of 



masses 



energetic particles @. 

Let us exemplify this simplified resummation method by recalculating the next-to- 
leading order screening mass in hot scalar electrodynamics. Separating static from non- 
static modes in i7oo(0, k) gives 

A,o(0,g) = e 2 ^{4ir 2 AA--2A} 

d 3 - 2e k 2 



e z Ta 



2= 



(2tt) 3 - 2£ k 2 + /i 2 

+ 2e 2 Ta 2 ^ £ J JL^JL {_ 4 [(* - q) 2 + ^]A n A n " + 6A n } . (4.14) 



n=l ' 



Here we have employed dimensional regularization in order to render the splitted expres- 
sions well-defined. From (|4.14 ) it is apparent that indeed only the n = contributions 
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are capable of producing a relative order e through e 2 T / fi — the other, nonstatic contri- 
butions can be expanded out in powers of fi 2 and will not give something nonanalytic in 
/i 2 , neither i.e., in e 2 . 

The next-to-leading order term is thus contained in 

d 3 ~ 2£ k 1 



W7 00 (0,g) 



2e 2 Ta 2e J 
2e 2 T> 



(2tt) 3 - 2£ k 2 + /i 2 



.,.„ T{-\+e) f4ira 2 ^ £ 



(4vr) 3 /2 I ^2 



+ 0(eV lnT ) 



O(eVlnT) 



e 2 7> 
2tt 



+ 0(ee^TfM) +0(eVlnT) , 



(4.15) 



where <r is the mass scale introduced by dimensional regularization. The limit e — > is 
regular so that there arise no UV singularities proportional to T, and the result coincides 



with (4.8). Actually, in the evaluation of 511, the use of dimensional regularization could 
have been avoided by subtracting off the unresummed zero-mode contribution, i.e. 

d 3 k 



5n 00 (0,q) = 2e 2 T J 



(27T) 



k 2 + fi 2 



1 

k 2 



(4.16) 



Compared to the full calculation in the preceding subsection, the g-independence of 
the next-to- leading order result is now manifest, which greatly simplifies its evaluation. 
The simplifications are less conspicuous in the case of 577^(0, q), which can be computed 
in an analogous manner, readily reproducing ( |4.10 ). 



5. Plasma frequency at next-to-leading order 



In this section we shall restrict ourselves to another limiting case, the one of long 
wavelengths, q — > 0, which again simplifies all calculations considerably. 

In this limit degenerates and contains only one independent structure function, 
because n 00 (Q ,q — > 0) = due to transversality and IIij(Qo,q — > 0) oc 5^ if is to 
remain regular. Since S^A^ = 5ijB l \ this entails that the longitudinal and the transverse 
polarization function coincide, II t (Q ,0) = IIe(Qo,0). In other words, without a wave 
vector there is no way to tell longitudinal photonic quasi-particles from transverse ones. 
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At q = 0, the unimproved one-loop result for T 3> Q ~ m reads 

n t ,e(Qo, 0) = ^ - ^»Q„ - In £ + 0(e 2 Q 2 T°) = m 2 - £-*Q m + 0(e 2 m 2 ) . 

(5.1) 

The mass m, often identified as the plasmon mass (although there is no longer any mo- 
mentum independent notion of mass) determines the plasma frequency, below which the 
medium cannot sustain free oscillations. At relative order e, ( |5.1|) implies a non-zero 
damping constant 

1 e 2 T e 

7 = -7^— Im n tAQo = m, 0) = — = — m . (5.2) 
2m 24n 87T 

However, we shall presently show that the subleading term is not stable under resumma- 
tion. In particular, the corrected 7 will turn out to vanish at relative order e. 



5. A RESUMMATION 

Proceeding as in the previous section we find for the resummed polarization functions 

4 



n t (Q o ,0) = i7 £ (Q ,0)=e 2 ^{--A; 2 AA--2A| 



r^ "(^> ( 2+ i /-Wf>M . (5 .3) 



3n 2 Jo VWTJF \ k 2 + fi 2 - {Q /2)\ 

Because of q = 0, the angular integration is trivial this time, and after separating off the 
constant T 2 -contribution, the 0(T)-part is easily evaluated, yielding 

n t (Q , 0) = n t (Q , 0) = ^f- + ^ j-/, - (y _ ( g / 2 ) 2 ]l - j . (5.4) 

At Qo — m this determines the correction to the plasma frequency, 

p 2 T 2 / S-v/2 1 - 9 \ p 2 T 2 
m 2 + 5m 2 = £i_ / 1 _ % j « £i_(i - .37e) , (5.5) 

as well as the 0(e 2 T)-contribution to the plasmon damping in the long- wavelength limit. 
However, (|5.3| ) is real at Qo = m, so that the damping constant vanishes at this order of 
magnitude. The bare result ( |5.2j ) is wrong because by not taking into account the thermal 
masses of the scalar particles a plasmon seems capable of decaying into those. In view 
of the long story of the plasmon damping puzzle in hot QCD H let us make the trivial 
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remark that the manifest gauge independence of the bare result ( |5.2| ) did not prevent it 
from being incomplete. The nonvanishing result for the resummed QCD damping constant 
by the way arises from the possibility of Landau damping, which is absent for thermal 
scalars. Thus in scalar QED an imaginary part to the polarization tensor appears only 
when pair decay becomes possible. Indeed, for Q > 2/i ( |5.4|) does become complex, but 
m < 2/i. 

In pure QCD at high temperature the next-to- leading contribution to the plasma 
frequency has recently been calculated by one of the present authors in Ref. |]IU| |, yielding 



(<5m 2 /m 2 )QCD ~ —Q.\%^/g 2 N'. Let us see how far this latter result can be understood 
by the above result on scalar QED. From the leading order terms it is clear that e 2 



corresponds to g 2 N, so we might try to apply Q5.4| ) by inserting the plasmon mass in 
place of the thermal mass of the scalars. This would give (5m 2 /m 2 ) « —0.028V g 2 N' , 
which is over a factor of 6 short of the actual result. Hence, the correction to the QCD 
plasma frequency is much larger than what might be expected from just the appearance 
of thermal masses in the loop integrals. 



5. B STATIC RESUMMATION 

In the preceding chapter we have seen that the next-to-leading order results on static 
Green's functions could be obtained also in a simplified scheme that resums only static 
modes. In the imaginary-time formalism, nonstatic modes are automatically hard so that 
their hard-thermal-loop corrections can be treated as perturbations. This scheme is par- 
ticularly advantageous in gauge theories, where it simplifies tremendously the calculation 
of static quantities and also of Green's functions with hard external frequencies. However, 
dynamical Green's functions with soft external frequencies are quite intractable because 
they require analytic continuation from the imaginary Matsubara frequencies which are 
either zero or hard. The separation into static and nonstatic contributions renders this 
analytic continuation to general soft external frequencies rather impossible. 

One might perhaps nevertheless expect that also in the nonstatic situation it is the 
resummation of internal static modes that gives the next-to-leading order correction which 
is nonanalytic in the coupling constant, and perform the analytic continuation on just their 
contribution. This is wrong, as we shall show in the present example of n t ,i(Qo, 0). 

Taking only the n = 0-contribution of the sum in (|5.3| ) and continuing at once to soft 
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Qo 7^ would give 



5n t>e (Q ,Q) 



-e 2 T 



4k 2 



d 3 k f 

(2vr) 3 \ (k 2 + /i 2 )(fc 2 + n 2 - Ql) ~ k 2 + /i 2 



(5.6) 



where dimensional regularization is understood (see above). This is readily evaluated, 
yielding 



5n t>i (Q ,0) 



e 2 T 



static contr. 



1- 








.Qo 





(Qg-/i a )+/i 



(5.7) 



It obviously disagrees with ( p.4[) . It has a wrong analytic structure in that the imaginary 
part sets in already at Qo > fi rather than 2/i, and of course does not reproduce the 
actual real contribution either. For instance, it would predict the correction to the plasma 
frequency squared to be 



m + 8m' 



static contr. 



e 2 T 2 ( _ 5yjT-9 > 



e 2 T 2 



9 



1 -0.09e) 



(5.8) 



which underestimates the next-to-leading order term by more than a factor of 4. 



Evidently, the nonstatic modes may not be neglected and even give the largest con- 
tribution in this example. (In Ref. ||25|| , the importance of the nonstatic modes was also 
noticed in the case of the gluonic plasmon damping.) The failure of the above reasoning 
may be traced to the premature analytic continuation (a collection of similar pitfalls with 



analytic continuation at finite temperature can be found in [2"6[|). In our case this analytic 
continuation was not possible because the polarization function in the imaginary-time 
formalism was given only at one point in the complex plane of soft energies, Qo = 0, 
which of course cannot be continued unambiguously into a function over nonzero (soft) 
frequencies. 



5. C CLASSICAL PHYSICS 

We close this chapter by trying to understand the results on the plasma frequency in 
the familiar intuitive terms of classical physics rather than through full-fledged quantum 
(field) theory. Looking at the hot scalar plasma as a simple system of particles interacting 
through Coulomb forces explains fully the leading order result u 2 = e 2 T 2 /9, as we shall 
see, and also part of the next-to- leading order result. 
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First, to count the states of scalar particles, imagine a box of volume V inside the hot 
plasma and apply periodic boundary conditions. Then 



N= (2^3 Jd 3 kn^p 2 + k 2 '^j (5.9) 

is the number of positive scalar particles in V as well as that of the negative ones. We shall 
consider again massless particles; the mass p in (|5.9Q is to make room for a dynamically 
aquired self-energy. Admittedly, the Bose distribution function n stems from quantum 
physics, but this is the only exterior element we shall employ. 

If a longitudinal electric plane wave E = E ~eisin(qx — cut) with an infinitesimal 
amplitude E is somehow activated in the medium, then Maxwell's equations tell us that 
there is no magnetic field associated with it. They reduce to 

E Q qcos(qx — out) = p , j\ = E u cos(gx — ut) . (5.10) 

In the long-wavelength limit q —>■ the charge density p vanishes. For the first component 
ji of the current density we may write j\ = 2e(N/V)vi, the factor of 2 coming from the 
oppositely charged particles (— e) moving with the opposite velocity —v±. By whatever 
Newtonian dynamics the Maxwell equations (|5.10|) are accompanied, it will lead to a 
factor A of proportionality between d t v and the force to the particles. Combining this, 
we obtain the plasma frequency u: 

d t v = XeE u 2 = X2e 2 ^ . (5.11) 



Guessing the inverse mass A to be simply 1/T, and determining N/V from ( |5.9|) as 
~ T 3 /8.2 leads to uo 2 m e 2 T 2 /4.1 — which is wrong. The simple resolution to this puzzle 
is that one has to take the relativistic version d t k = eE of Newtons equation and identify 
the velocity v with a mean value: 

where / is the distribution function (which in equilibrium is the Bose function) and 
ki I \J p, 2 + k' 2 ' is nothing but the relativistic velocity- momentum relation. Particles when 
accelerated carry their probability with them, so 



d t f = -V*/ d t k = -eE x d kl n dp 2 + k 2 + 0(E ) . (5.13) 
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Using this in (|5.12|) and ( |5.11|) we end up with 

IV 1 r „ k 2 



which when taken at /x = gives A = (V/N)T 2 /18. Thus we obtain indeed cu 2 = e 2 T 2 /9, 
where one factor 1/3 came from averaging over the directions of k and the other from 
the integration in (|5.14|) . 

Trying now to go beyond leading order, we could take into account the thermal mass 
fi = eT/2 acquired by the scalar particles. From formula (|5.14 ) we then obtain for the 
'classical' next-to-leading-order plasma frequency 

2 r - U (yk T T~jj?' ) I 



Classical " j dk k " JWTW ~ ( 2 + ~k^T7 ) ' (5 ' 15) 



which is to be compared with the true result (|5.3| ). The only difference is the missing 
back-reaction in the form of corrections involving Q = o; c i ass i ca i in the integrand. In ( |5.3| ) 
these terms give rise to an imaginary part when Qq > 2/i, which obviously corresponds 
to pair creation. This cannot be captured in purely classical terms, of course. 

Despite the missing terms, ( |5.15|) is extremely close to the correct result (|5.5| ), to wit 

e 2 T 2 / 9 \ e 2 T 2 
m 2 + 5m 2 | c i assica i = — — 1 - — e w — — (1 - 0.36e) . (5.16) 

9 V 87T / y 



6. The complete plasmon spectrum at next-to-leading order 



Because of the simplicity of the resummed scalar propagator determining the two 
polarisation functions Hi and II t up to and including the next-to-leading order terms, 
these can in fact be calculated analytically for the entire cu-g-plane. In contrast to the 
QCD counterparts, the spectral density of the dressed scalar propagator A = 1 / (K 2 — fi 2 ) 
has no cut contribution: 

r l l 

A(K) = dx— p(x,k) , p(x,k) = — [A(x-u) - A(x + aJ) , (6.1) 

J J\r\ — x 2lo 
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where u = \J p 1 + k' 2 ' . II \ and lJ t as given in (|3.5| ), fl3.6|) and 
two sums (and integrals) 



and 



> g) 2 



are determined by the 



(6.2) 



Here, the upper index on A refers to the shift Kq — > Qo — Kq as well as to k — > g — . 
Hence, the propagator A - has the spectral representation ( |6.1| ) with K shifted and with 
the frequency cU_ = \j \j? + ( k — q) 2 ' in place of cJ. Equivalently, we may view it as 
'another propagator' taken at Q — K and with a spectral density p_ as described. At 
this point we recall two formulae derived earlier, (6.6) and (6.8) in Ref. [ 10| . They concern 
the real and imaginary parts, respectively, of expressions like (|6.2|) . The two are related 
by the dispersion relation 



Ke^A"A/(^; 



dt 



1 r 



SmJ2 A~A/(fc 



UJ=t 



(6.3) 



We may therefore concentrate on 
SmV A"A/(jfe) = -kuT (— ) fd 3 kf(k) [dx—^ -p(x, k)p(x-uj, k) , (6.4) 

V27T/ J J x{x — UJ) 

where on the right-hand-side only the leading temperature-dependent part resulting from 
n(x) ~ T/x has been included, assuming that the integrations are restricted to soft 
arguments of the Bose function either automatically or after suitable subtractions. 



The details of first evaluating the simpler imaginary part (|6.4|) and then using ( |6.3|) are 
given in the Appendix B. The results may be put together as follows. The prefix 5 again 
indicates that the known leading hard-thermal-loop contributions have been subtracted. 
Three regions in the u;-g-plane are to be distinguished: uj 2 < q 2 (region I), q 2 < u 2 < 
4/j, 2 + q 2 (region II), and 4/z 2 + q 2 < uj 2 (region III). Then: 



5TI e 



Silt 



e 2 T u 2 -q 2 ( . _ u 2 



■in 



4/i + 2z£ [TZ + iJ 

Q 



e 2 T 



•57T 



-8fi + 



4/i 2 + q 2 



to 



q 



[K + iJ] 



(6.5) 
(6.6) 



(6n g - 6n e ) 



2 
16tt 



-4/i 



uj 2 + q 2 



to 



-2iS + 



Afi 2 



(6.7) 
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Figure 3: The real part of bUi in units of e(eT) 2 . LU and q are in units of eT. 



where 



K 



(Q. + q\ fQ-q" 

arctan — arctan 



2/i 



2/i 



2 arctan 



Q 



J 



£ 




2/i+ \Q 
ttq + uo 2 



in I and III 
in II 



fig — uo 2 



in I and III 
in II 



in I and III 
in II 



(6.8) 

(6.9) 
(6.10) 



with fi, if real, the positive square root of 

uj 2 - q 2 - 4/i 2 



n 2 = u 2 



uj 2 — q 2 



(6.11) 



As one learns in Appendix B, the combination IZ + iJ can be cast into the compact form 

2/i — i£ — iq\ 



Tl + iJ = i In 



2/i — i£ + iq 



(6.12) 



The static limit of these results obviously reproduces the ones derived before, 
and ( [4.10 ), whereas a bit of calculation is required to verify that the long- wavelength limit 
q — > indeed gives SIIi(Q , 0) = SII t (Q , 0) and coincides with the result derived in (|5.3|). 



Notice that 5IIi and SII t are purely real in region II. In region I, where the hard- 
thermal-loop contribution has an imaginary part corresponding to Landau damping, there 
is now a correction term of relative order e, whereas the imaginary part appearing in region 
III is the leading term resulting from the possibility of the decay of virtual photonic plasma 
excitations into pairs of scalar quasi-particles. 

The real part of 5IIi and 5II t is displayed in figs. 3 and 4, respectively. For a com- 
parison to the leading-order result see Fig. 2. In Sect. 5 we have seen that the effect 
of the next-to-leading order contributions is to lower the plasma frequency according to 
(p.5|). We shall now study the whole spectrum of propagating photonic quasi-particles, 
i.e. g / 0. 
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Figure 4: The real part of bTLt in the same units as used in Fig. 3. 



In the case of transverse photonic excitations, it turns out that for increasing q the 
ratio of the corrected frequency uj t {q) to the lowest-order one decreases somewhat. For 
q ^> eT the effective thermal mass becomes momentum- independent, and the dispersion 
curve goes over into a perfect mass hyperboloid with asymptotic mass 

< + Sm ~ = ^f(}-^) ■ ( 6 - 13 ) 

We now turn to the the more interesting longitudinal plasmons, which are collective 
modes without analogues at zero temperature. As we have seen in chapter 2, their leading- 
order dispersion curve approaches the light-cone exponentially with increasing q. At the 
same time the residue of the corresponding pole in the propagator vanishes exponentially 



271, quickly rendering them unimportant for larger values of q. 



Including the next-to-leading order terms now, we expand everything around the 
leading-order result u)o(q) 2 and rewrite the condition u 2 = q 2 + II^u, q) for fixed q as 



u 2 = oo 2 + 5II e (u , q) + (u 2 - lo 2 ) c^iirVo, q) + 0(e 2 m 2 ) (6.14) 

which gives 



6u = u -u = - a - nhard , : ^ + 0{e m ) . (6.15) 



5II e (uj , q) 

In order to work consistently at relative order e, we drop all contributions that are of 
higher order. With Z7] iard given by ( |2.6|) , ( |6.15|) can be rewritten as 



2uj 2 

cu 2 = q 2 + i7, hard K, q) + - 2-L ° 2 5 6II e (u , q) . (6.16) 

The resulting dispersion curves for e = 0.3, 1, and 2 are shown in Fig. 5. It is seen 
that the corrected curves are not just a slightly down-scaled version of the leading-order 
result. The dispersion curves now hit the light-cone at a finite value of q, above which 
there is no longer any solution corresponding to propagating plasmons, since for uj < q, 
Hi has a large imaginary part ~ e 2 T 2 , which prevents the existence of weakly damped 
excitations. 

The residue of the plasmon pole can be defined in a gauge-independent manner by 
projecting the propagator G flll (Q) onto conserved currents J with Q^J^ = and restricting 
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Figure 5: Leading and next-to-leading order dispersion curves of the longitudinal plasmons. The upper 
curve is the leading-order result; below it are the corrected ones for e = 0.3, 1, and 2, respectively. The 
version on the left consistently discards contributions to 5(jJ 2 (q) that are beyond relative order e; the one 
on the right gives the location of the poles in the propagator when all next-to-leading order contributions 
to the self-energy are kept. 

Figure 6: The residues Zf(q) of the plasmon poles associated with the dispersion curves of Fig. 5. 

to spatially longitudinal J} = P . From the component of the propagator associated 
with Jj in 



Q 2 -?* 

we extract the residue as 



J fl G flu (Q)r = -^J'A e + ... , (6.17) 



Z t = lim . 6.18 

Other definitions, differing in normalization, are possible [J27J], but the one given here 
avoids unnecessary kinematical singularities. Abbreviating $ = UqIIi(u , q) / (u 2 — q 2 ), 
( |6.18|) gives 

z t ( q ) = z'r d [ i + z^ d (a^m + 5u 2 di^r d ) } , (6.19) 

where Z} &rd = 1 - 9 w 2$ hard . 

In Fig. 6 the leading-order result and the corrected one are compared for the same 
values of e as in Fig. 5. For both, the residue decreases rapidly for increasing q, with 
the one-loop-resummed result being smaller than the leading-order one for all q > and 
arriving at zero at a finite value of q. 

Both results indicate that there is now only a finite range of q and lu in which lon- 
gitudinal plasmons can exist. However, the very fact that the corrected result changes 
qualitatively signals that a break-down of the resummed perturbation theory is occurring. 
This qualitative change is possible only by 5IIi/Q 2 becoming greater than Ilf^/Q 2 , 
although the former is down by one power of e. The reason is that n} wA /Q 2 is logarith- 
mically divergent for Q 2 — > 0, but 5IIg/Q 2 diverges even stronger, like 1/y/Q' 2 ' (see ( |6.5| ) 
and (|6.1CI|) ). Hence, y/Q' 2 ' /q ~ e eventually (over-) compensates for the factor of e in 5IIi. 
In fact, had we solved Q 2 = i7^ ard + 511 ^ for a small finite value of e without expanding 
ui{q) = u>o(q) + eoj\{q) + ... as we have done above, the result, shown in the right half of 
Fig. 5, would have been that the dispersion curve follows the result uio(q) + eui(q) closely 
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until Q 2 /q 2 becomes comparable to e, after which it bends back sharply, running down 
hard by the light-cone back to u = q = 0. If the higher-order corrections to IIi/Q 2 keep 
being more and more singular for Q 2 — > 0, this kinky behavior may change from order to 
order. It is clear that at any finite order one could trust the result only up to a certain 
small distance from the light-cone. 

In fact, this break-down of our resummed perturbation theory is not actually due to 
higher-order diagrams, but is caused by the break-down of the high-temperature expansion 
of n^ u at Q 2 = 0. Since n e /Q 2 = —II 00 /q 2 and the internal lines in 77^ have been dressed 
by their thermal masses, a kinematical singularity at Q 2 = should not be there at all. 
Indeed, as shown in Appendix C, performing the high-temperature expansion directly at 
the light-cone Q 2 = 0, i.e. evaluating 77 Mi ,(Qo — <?;<?)> the leading and next-to-leading 
terms are found to be 



iTfsum. e 2 T 2 



lim / v , , „ 

Qo^g Q 2 3q 2 



2T 1 C'(2) 



In h - - 7s + 



/x 2 '* C(2) 



+ O(eVT°) (6.20) 



with being Euler's constant and ( the Riemann zeta function. The first term on the 
right-hand- side is ~ ln(l/e), so the original logarithmic singularity evidently got cut off 
at \JQ 1 ' /(eT) ~ e by the resummation of /i ~ eT. The equation 1 = FIi/Q 2 with Q 2 — ► 
therefore has the solution 



^./(eT) 2 = ~ 



, 4 1 C'(2) 



e 1, 2.094... e 
= - In . 6.21 

4tt 3 e 4tt v ' 



Because higher-order diagrams are not singular at Q 2 = 0, either, the result (|6.21 ) is stable 



up to and including order e. The fact that it is non-analytic in e makes it clear that it 
could not be obtained by organizing the resummed perturbation theory in powers of e. On 
the other hand, the strictly perturbative next-to- leading order result, depicted in the left 
half of Fig. 5, becomes inaccurate only close to the light-cone. Indeed, the "perturbative" 
result for g crit w 0.77eT for e = 0.3 is not far off from the one of Eq. ( p.21| ), which gives 
9cnt. ~ 0.79eT. The alternative result in the right half of Fig. 5, however, where 5u 2 was 
not truncated at relative order e, but the entire next-to-leading order result for 5IIe was 
kept to define a new full propagator, turns out to be qualitatively wrong. 

The fact that the resummed result for 77 \j Q 2 is no longer singular at the light-cone also 
implies that the residue of the plasmons does not completely vanish there. This means 
that there are longitudinal massless photonic excitations with a fixed value of q = uj. 
For larger values of uj and q, Q 2 becomes negative, i.e. space-like, and Landau damping 
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sets in. At the level of hard thermal loops, the imaginary part ~ e 2 T 2 switches on 
discontinuously, Qmlli^Q 2 ) / Q 2 = 9(— Q 2 )^-m 2 uj/q 3 . Because the logarithmic singularity 
in the hard thermal loop is removed, this discontinuity is smoothed out. As we show in 
Appendix C, for to < q and q 2 —uo 2 <ti e 2 g 2 , an extra factor exp(— e^J 8 ^ q 9 _^ ') arises — the 
imaginary part sets in with all of its derivatives vanishing. There is now a finite range in 
the space-like region with small damping. Consequently, the plasmons are removed from 
the spectrum only for (q — u)/(eT) > e 2 . 

This phenomenon of a finite range of q and u, where longitudinal plasmons can exist, 



was predicted in the case of QCD in ER B9f. The simplicity of scalar electrodynamics 



has allowed us to study it in full detail. In the case of non-ultrarelativistic QED (i.e. 
with electron mass m > T), a finite range for the longitudinal plasmons has been found 



already in 1961 by Tsytovich [30]. 



The simple step beyond the resummed perturbation that was necessary in the above 
calculation of q cvl t. might perhaps shed some light on how similar failures with Braaten- 
Pisarski resummation could be overcome. In Eq. ( |6.20| ) we have witnessed the hard 
thermal loop itself becoming subject to change under resummation of the hard thermal 
loops close to its singularities. So there are cases where the resummed perturbation theory 
cannot be organized in excess powers of e. A perhaps similar failure was encountered in 
a recent attempt to calculate the production rate of real non-thermal photons in a QCD 
plasma by Braaten-Pisarski resummation . 

We close this chapter by suggesting that such kinematical singularities in the Braaten- 
Pisarski scheme for QCD can be taken care of by extending resummation to the hard 
thermal loops themselves in the same way as we have done above for scalar electrody- 
namics. At hard loop momenta it is in fact not necessary to use the full complicated form 
of resummed propagators and vertices in QCD for deriving the modification of the hard 
thermal loops near the kinematical singularities they otherwise give rise to. The latter are 
caused by the masslessness of the propagators in the bare hard thermal loop, so it will be 
the true quasi-particle poles (rather than cut contributions) that are important. At large 
loop momenta 3> gT, and only there, these are given by momentum-independent masses 
(cp. Eq. ( |6.13| )), whereas the additional collective modes have exponentially vanishing 
residues at leading order. Only at relative order 0(g) the form of the dressed propagators 
and vertices at soft momenta become relevant, whereas the leading terms are determined 



by simple massive loop integrals. The QCD result analogous to the part oc T in ( 6.20 ) 



has the same form for the bosonic (gluon) contributions with e 2 — > g 2 N, fi — > moo. In 
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the case of fermionic (quark) contributions a similar formula holds, obtained by replacing 
7_b - > 1e — In 2 and taking now the asymptotic value of the thermal quark mass. With 
Nf flavors, and iV + = N + Nf/2, this yields at leading order 



?crit. 



3m 



N 2^ 



+ 



2N: 



In 



8 1 _ C(2) 
gy/C? 2 lE+ C(2) 



(6.22) 



where (7p = (A^ 2 — 1)/(2N). (This result differs somewhat from the one of Ref. [28|, but 
this seems to be due simply to an incorrect evaluation of the integrals given therein, with 
which we otherwise agree.) We intend to cover the case of QCD more fully in a future 
publication. Let us just note here that g 2 rit as given by (|6.22|) goes down to zero when 
g = 1.48 ... (for = 3, N f = 0) or g = 1.31 ... (for = 3, N f = 2). Clearly, next-to- 
leading results are becoming decisive here to determine the fate of the QCD plasmons. 



7. Resummation of scalar damping 



Up to now we have studied exclusively the effects of the thermal mass fi acquired by 
the scalar fields, i.e. of the leading term of the self-energy of the scalars. A calculation of 
the next-to-leading order terms in this self-energy would already involve the more compli- 
cated thermal photon propagator, making it necessary to resort to numerical integrations. 
However, the limiting case of large scalar momenta turns out to be accessible by purely 
analytical means. 

For external scalar momentum Q with Q ~ T, the next-to-leading order terms can 
be most easily derived by using the static ring resummation (see Sect. 4B). For this we 
only need to resum the zero-mode propagators 

,(0,*) = ^ , G (0,*) = ^ + (a -l)^-_^ (7 ,) 

of scalar and gauge fields, respectively, with the four- vector K = K — (KU)U and = 
e 2 T 2 /3, see (^9|) , (|2~T0|) and Q. 



The correction term to the scalar self-energy, for large Qo, q ~ T, is then given by 
Sa = j2^J dk ^G, V {K) ^—y 2 G, V {K) 



25 



K 2 







(7.2) 



J K =0 



where, while the first term clearly corresponds to the tadpole diagram, the loop contributes 
twice, once with the hard momentum Q running through the scalar line (second term) 
and once with Q in the photon line (third term). Only the propagator that carries the soft 
momentum K requires resummation. On mass-shell, Q 2 = fi 2 , the real part of ( |7.2|) gives 
the next-to-leading order term in the thermal mass of energetic scalar particles, yielding 
(after some tedious integrations) 



fy 2 = — sfte 55 



e 2 T ( 1 
m eL + —fJ, 



2tt 



1 



2tt Iv^ 



1 



(7.3) 



The imaginary part, which determines the damping rate ( of energetic scalar particles 
defined by 



55 



2QoC 



(7.4) 



however leads to an infrared-divergent expression, which is in fact virtually identical to 
what is obtained in QED or QCD || 
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-(1 -a) 
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+ 0(A L 



(7.5) 



There is a logarithmic singularity caused by the absence of screening for static transverse 
electromagnetic fields. In non-Abelian theories the infrared singularity in (|7.5| ) could 
perhaps be removed by the appearance of a magnetic screening mass, but in the Abelian 
case i7n(0, k) oc k 2 . In Ref. |29|, [32|, Q it was suggested that this infrared divergence is 
instead cut off by the damping of the energetic particles itself, entailing A ~ e 2 T and 



c 



e 2 T 1 
In - 

47r e 



(7.6) 



However, this cut-off depends on having the external momentum slightly off the assumed 
complex pole |34J and this in turn gives rise to gauge dependences at the order e 2 TA°. 
There is still an on-going discussion on how to remedy this situation and to become able 



to go beyond the leading logarithmic term of ( |7.6|) , see |p5 |. In particular it seems likely 
that at least in Abelian theories the propagators of moving quasi-particles do not have 
simple poles on the unphysical sheet but rather branch singularities. 
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At any rate, in order to go beyond the approximations leading to ( |7.5| ), one evidently 
should include self-consistently physics at scales below eT. Taking the result of (|7~"6| ) for 
granted, one has ( ~ /ieln(l/e) ^> 5fi, so seemingly damping effects are the next impor- 
tant corrections to be included in an improved resummation scheme. The consequences 
of this have been analysed in Ref. []32| in the case of QED, and recently also in scalar 



electrodynamics in Ref. | 12| . Here we shall follow the lines of Ref. |32| and correct what 



has been presented in Ref. [12]. 



If one tries to improve the propagator for scalars at large momenta by including the 
unexpectedly large damping (|7.6|) , one finds that in contrast to the thermal mass a small 
but finite damping constant may even change the leading hard thermal loop result for 



small momenta. It does so, however, by violating gauge invariance [j3~2~ |. In particular one 
obtains 

2iC e 2 T 2 

which is inconsistent with transversality of the self-energy. This seems to have gone 
unnoticed in a recent paper on the analytic properties of finite-temperature scalar elec- 
trodynamics Ref. P0|], where it was claimed that inculsion of damping would restore 



analyticity of the zero-momentum limit of thermal Green's functions. 

In order to include a finite damping self-consistently one obviously has to consider 
vertex corrections. Lacking an effective action, which would furnish corrected Feynman 
rules and compensating ("thermal") counterterms, one may analyse the possible contri- 
butions from vertex functions through the coupled Schwinger-Dyson equations they have 
to satisfy. The Schwinger-Dyson equation for the photon self-energy is depicted in Fig. 
|7|. It involves the fully dressed scalar propagator as well as bare and dressed vertices. In 



Ref. [32, 12] it has been argued that in the infrared the dominant corrections coming from 
the vertices are given by ladder diagrams and that their leading terms can be obtained by 
solving the Ward identities of vertex diagrams in terms of the self-energies and discarding 
possible transverse contributions. 

In order to avoid the intricacies of analytic continuation, in particular that of vertex 
diagrams, we shall like Ref. P2L |T!J work in the real-time (Schwinger-Keldysh) formalism 



36(1 • hi this formalism perturbation theory is based on a time contour that runs along 
the real axis and back, which leads to a doubling of the fields corresponding to the part 
of the contour they are living on, and to a contour-ordered propagator that is a matrix 
with respect to this bisection. If a, b = 1,2 denotes the field type, a scalar propagator is 
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Figure 7: Schwinger-Dyson equation for the photon self-energy. The lines with a black bullet are the 
resummcd propagators of figure 1. 



now given by the matrix 

- iA(x -y) = (T c <j) a (x)<j) b (y)) 



(T<j>(x)<i>(y)) (0(2/)0(x)> 
(</>{x)<f>{y)) (T<f>(x)<f>{y)) 



{7.1 



where the contour-ordering T c is broken up into time-ordering T for type-1 fields and 
anti-time-ordering T for type-2. Mixed propagators correspond always to a fixed order of 
the operators because type-2 fields are always further down the contour than type-1. 

Analytic continuation from the imaginary-time formalism eventually leads to a col- 
lection of Green's functions which depend on the prescriptions chosen for each external 
line. In the real-time formalism, these are given by the various linear combinations of the 
components of the corresponding matrix quantities [37]]. The retarded and advanced 
propagators are given by 



A R (P) = A n (P) - A 12 (P) , A A (P) = A n (P) - A 21 (P) . 
The symmetric combination, defined by 

A P (P) = A n (P) + A 22 (P) = A 12 (P) + A 21 (P) , 
is, in thermal equilibrium, related to the former, 

A P (P) = (l + 2n(P ))(A R (P)-A A (P)) . 



(7.9) 



(7.10) 



:7.U) 



Whereas at tree level, vertices are exclusively connecting either type-1 or type-2 fields, 
there is a large variety of full vertex functions which differ in their analytic properties. In 
the first diagram on the right- hand- side of Fig. [7| we have a full vertex function connecting 
an incoming scalar, a photon, and an outgoing scalar. Defining retarded vertices with 
respect to these three lines in turn we have 



R(i) 



a.b 



R(7) 



a.b 



R(o) 



a.b 



(7.12) 
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where the indices are in the order of incoming scalar (i), photon, and outgoing scalar (o). 
There are various symmetric combinations, of which we shall need only 



r P = E r ^ ( 7 - 13 ) 

0,6 

which is the one that singles out the photon line. 

With these definitions the retarded piece of the first diagram on the right-hand-side 
of Fig. [F] is given by 



' 1 K -{2K + Q),{A R (Q + K)T^[A P (K) + A R (K)] 



2 J (2tt) 4 

+ [A P (Q + K) + A A (Q + K)]T^A A (K) 

+ A R (Q + K)T P U A A (K)} . (7.14) 

The second diagram involves only the bare 4-vertex and therefore differs from the 
usual tadpole (or sea-gull) diagram only in that the full propagator has to be taken. In 
terms of the real-time quantities it is simply given by 

n$ R (Q)=ie\ v J 0^[A P (K) + A R (K) + A A (K)} . (7.15) 

The other diagrams are of higher (explicit) power in e and we shall not need them in the 
following discussion. 

We wish to consider a resummation of both thermal masses and the "anomalously" 
32| large damping of the scalars. We therefore begin by adding the damping to the 



retarded and advanced scalar self-energies according to 

- E R {P) = i?- 2tP ( , -S A (P) = /i 2 + 2tP ( = -~* = -S R {-P) . (7.16) 

The momentum dependence of these additions spoil the Abelian Ward identities and this 
leads to a violation, proportional to (, of the generally valid transversality of the photon 
self-energy already at the order of e 2 T 2 , i.e. at the level of hard thermal loops [32|. In 



Ref. |32|] it has been shown however that higher-order vertex diagrams start to contribute 
because with the above modification of the propagators there are now contributions pro- 
portional to £ from would-be pinch singularities which are now cut off by the small but 
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finite damping. It turned out that the dominant vertex contributions just correspond to 
solving the Ward identities according to 



-ie{2K + Q) l + 



T p (K + Q,K) 



-ie(2K + Q). 



ie(2K + Q), 



1 + 



2KQ + Q 2 

giCgo 
2KQ + Q 2 

1 



[E R {K + Q)-E R {K)[ 



(7.17) 



E P {K + Q)- E P (K)\ 



,fl 2KQ + Q 2 
-ie(2K + Q) M r 

~2kq + q 2 ~ { C) r° 

+ 2(K + Q )n(# + Qo) - 2K n(K ) 



(7.18) 



(In the case of QED [|32| only T p needed correction because there the damping of the 
electrons arises from a constant contribution to the self-energy.) 



Inserting the corrected vertices into ([7.14 ) and keeping only terms involving the distri- 
bution function n, one finds for the dressed one-loop contribution to the retarded photon 
self-energy 

d A K 



n*{Q) = -2ie 2 j —{2K + Q),{2K + Q) U { 



A R {K + Q) 



1 + 



2<Qo 



2KQ + Q 2 



n( 



(K ) (A R (K) - A A (K)) 
2K Q + Q 2 **{K + Q)*(X)} 



+ 2< (Ko±Qo)n(K + Qo) - K n(K ) aR 



(7.19) 



By a change of variables, this can be brought into a form that can be more easily 



compared with the corresponding expressions given in Ref. [12 



n*(Q) = 2te 2 J^(2K + Q),(2K + Q) u n(K ){ 

4<K 



(27T) 

A R (K + Q)A R (K) 



2KQ + Q 2 
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A R {K + Q) (A R (K) - A A (K) 



l + 2< 



2K + Qq 
2KQ + Q 2 



+ ^e 2 g,„J ^n(Ko) (A R (K) - A A (K) 



(7.20) 



again up to temperature-independent contributions. This does not completely agree with 
Eq. (25) of [I2|: there the first term in the first integrand was omitted, and in the second 
integral a dressed 4-vertex was used, contrary to what the Schwinger- Dyson equation 



prescribes. Moreover, it was assumed in Ref. |12[ that the dominant contributions from 
the vertices would be the correction terms, because they involve a factor 1/(2KQ + Q 2 ), 
so the tree-level vertices were deleted. Further, (2K + Q) /X /(2KQ + Q 2 ) was simplified 
to K^/KQ. Neither of these truncations is justified, however, as a complete calculation 



readily shows. The part of ( |7. 14| ) which involves only the tree vertices does produce the 
same leading temperature results as without resummation of the scalar damping plus 
terms proportional to (, and the vertex corrections that are proportional to ( stay so 
when the integrals are evaluated. The replacement {2K + Q) fl /{2KQ + Q 2 ) — > K^/KQ 
on the other hand changes the contributions coming from the vertex corrections, but even 
after discarding the pure tree-level contributions the leading temperature contributions 
are not restored. For instance, the Debye mass iloo comes out minus half the correct 



value. In Ref. |]T2| , the missing factor of 2 was compensated by taking the first diagram 
in the right-hand-side of Fig. [5] twice; the minus sign seems to have been simply lost. 



Another thing that goes wrong with the expressions put forward in Ref. [12 is in fact 
gauge invariance. As a special case of the transversality of the photon self-energy one 
has to have II 00 (Q ,0) = 0. The expressions given in Ref. |12| violate this requirement 
already in the leading terms of the high-temperature expansion, both with and without 
the unjustified simplifications performed therein: II 00 (Q ,0) in each case turns out to be 
proportional to ie 2 T 2 (/Q . 



Coming back to the complete results ( |7.19|J7.20| ), they can be evaluated by closing the 
contour in the .fTo-integration. In the upper half plane, there are poles from the advanced 
Green fucntion G A (K) and an infinite set of poles from the Bose distribution function 
n. There are also potential poles from the denominators of the vertex corrections, which 
are ambiguous because the Ward identities do not determine the prescription for these. 
We shall therefore consider only such cases where these denominators cannot give rise to 
poles on algebraic grounds. This holds generally for Q = and with Q ^ at least 
for i7 o(Qo,0), which is of special interest as a non-zero result indicates a loss of gauge 
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invariance. 



In the form ([7.20 ) there are nonvanishing contributions from the poles of the Bose 
distribution function. For example, with ~q = 0, Qq <C T, they contribute at leading order 
in T 



- 16e 2 TC £ / 
ie 2 (T< 



n=l ■ 

,2/^2 



(2tt)( 3 - £ ) iQ ((27mT) 2 + P) 2 



3g c 



W + )+0(e) , (7.21) 



which is due to the first term in the first integral of ( 7.20Q . The omission of this contribu- 
tion is what led to the (unnoticed) loss of gauge invariance in Ref. \T2\. In the seemingly 
more complicated form ( [7.19|) , however, it turns out that in all the considered cases the 
contributions from the poles of the distribution functions which would be either propor- 
tional to T 2 or T exactly cancel upon integration over spatial momenta; only the poles of 
G A (K) contribute by closing the contour around the upper half plane. 



At leading order T 2 , the correction terms to the vertices in Eq. ( 7.19 ) are found 
to exactly cancel the result ( [7.7|) which indicated that dressing of only the propagators 
violates gauge invariance. However, contrary to what has been claimed both in Ref. p2 
and [12|, the one-loop resummed expressions are not completely independent of (. There 
are still contributions oc T in (|7.19 ) that are dependent on (. Checking gauge invariance 
first, we find for the leading term that violates transversality of the photon self-energy 

ie 2 C 3 T 

^oo(Qo,0) = ln \ or , ■ (7-22) 

vr/i(g + 2<) 

For Qo<C-> this is of order e 5 ln 3 (l/e)T 2 ; for larger Q$ it is of even higher order. (A 
similar result is obtained also in the case of QED.) 

Inspecting now potential contributions to 11^(0, q), we find that the only new correc- 
tions are of precisely the order where gauge invariance is being lost, 

ffM(0 ,^0) = £f-£^-^, (,23) 

n«(0,g-0) = ■ (7.24) 

The latter result, which if correct would correspond to a magnetic screening mass, is by 
the way just another expression of the fact that gauge invariance has been lost — in 
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scalar electrodynamics one can easily prove the absence of a magnetic screening mass to 
all orders of perturbation theory [[U], |38[1 . 

The form of the above results in fact shows why the resummation of scalar damping 
is bound to fail eventually. The contributions involving ( are obviously dominated by 
the infrared region of integration, which is effectively cut off by the thermal mass /i. The 
damping of the scalar fields on the other hand has its simple constant form certainly only 
for momenta 3> eT, so precisely where interesting effects might show up, this resummation 
is clearly insufficient. 

In view of recent attempts to go beyond the Braaten-Pisarski resummation scheme 
by a simple resummation of damping effects f39 |, we close by emphasizing again that a 



resummation of damping only in the propagators proves to be incorrect and insufficient 
already in the Abelian case. Vertex corrections are of vital importance to keep gauge 
invariance, and they tend to just undo the resummation of damping in the propagators. 
Moreover, in the infrared regime where a resummation of damping might eventually be- 
come important, the full-fledged self-energy and vertex functions will be relevant, which 
can hardly be approximated by a constant damping term. 



8. Summary and conclusions 



The simplicity of hot scalar electrodynamics has allowed us to determine the complete 
next-to-leading order corrections to the spectrum of photonic quasi-particles. For the 
most part, this was a straightforward application of the resummation scheme of Braaten 
and Pisarski developed for QCD, with the bonus of being able to do all calculations 
analytically, so that we did not depend on numerical approximations. As we have shown 
in detail in the Appendix, these calculations can be most efficiently done by evaluating 
first the simpler imaginary parts and exploiting dispersion relations to obtain the complete 
expressions. 

Against these results we tested the simplified resummation scheme of Ref. |J , which 
is a systematic extension of the old ring resummation prescription of dressing zero modes 
only, and found that it works only for static quantities (for which it has been put forward 



33 



originally), whereas non-static Green's functions turn out to receive next-to- leading-order 
corrections also from non-static modes. In the case of static Green's functions, on the 
other hand, the resummation scheme of Ref. constitutes a computational simplification, 
which is of course more pronounced in non-Abelian applications. 

The next-to-leading order results for the plasmon spectrum constituted mostly small 
corrections to the ones derived from a bare one-loop calculation, but on two occasions 
they led to a qualitative change. 

First, the bare one- loop result gave a nonvanishing result for plasmon damping of the 
order of e 2 T ~ em. Like in QCD, this result turned out to be inaccurate, although unlike 
the bare one- loop gluonic plasmon damping, it was positive and gauge- independent. The 
resummed result revealed that the photonic plasmon damping is in fact zero at relative 
order e, because the thermal scalars are too massive to be produced by pair creation and 
also because at this order there is no Landau damping contribution from scalar quasi- 
particles. The obvious moral is that gauge-independence is only a necessary criterion for 
a correct result, not a sufficient one. 

Second, the next-to-leading order results for the longitudinal plasmons turned out 
to break down very close to the light-cone. At leading order, the longitudinal plasmon 
branch of the dispersion curves approaches the light-cone exponentially. Very close to 
the light-cone, the next-to-leading order corrections derived strictly along the lines of 
the Braaten-Pisarski resummation scheme begin to dominate over the leading-order con- 
tributions, so this scheme ceases to be actually perturbative. The reason for this is a 
singularity, at the light-cone, in the hard thermal loops that are being resummed. A 
more persistent resummation removes this singularity and allowed us to calculate the 
true next-to- leading order dispersion curve, which hits the light-cone at a finite value 
of momentum ~ eTln 1//2 (l/e), and continues even a short distance below the light-cone 
before Landau damping prohibits the existence of weakly damped excitations. We have 



argued that the same phenomenon occurs in QCD, as has been done before in Ref. [28 
whose corresponding result we believe to have corrected. The encountered break-down of 
the canonical Braaten-Pisarski scheme and its redress could be relevant also with respect 



to another recently observed break-down of this scheme in Ref. ||31|| , which also is due to 
an irregular behavior of the hard thermal loops at the light-cone. 

Finally, we have investigated another modification of the Braaten-Pisarski scheme, 
which attempts the next step by resumming also damping contributions from internal 
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lines. In accord with Ref. [p2| , where similar issues have been studied in the case of 
fermionic QED, we found that there are no contributions to the photonic polarization 
tensor from such a further resummation, if and only if the vertices are corrected in addition 
to the propagators. We disagree with the findings of Ref. |12|], however, where this sort 
of resummation was examined for scalar electrodynamics. In contrast to Ref. []12|, we 



found that the corrections to the vertices do not dominate over tree contributions, an 



assumption which led Ref. [12] to postulate a new diagrammatic framework. Instead, both 
contribute at the same order of magnitude and give correct results with the conventional 
set of diagrams as prescribed by the Schwinger- Dyson equations. Further, we observed 
that the assumption of a constant damping term can be sufficient only for leading-order 
results, and that it leads to a violation of gauge invariance at higher orders even with 
vertices corrected such that they satisfy the Ward identities. However, without correcting 
the vertices, gauge invariance is lost already at lowest order. This seems to have gone 



unnoticed in a recent paper on analyticity properties of scalar electrodynamics, Ref. [[40 
It also casts doubt on recent attempts to improve on perturbative results by an ad hoc 
resummation of damping contributions in QCD f39[ only in propagators. 
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Appendix A 



In this Appendix, which supplements Sect. 3, we briefly go through the construction of 



the effective action summarizing the hard thermal loops along the lines of Refs. [^2, [II, [42 
restricted to the case of scalar electrodynamics. As was first pointed out in p2|, this 



effective action is completely determined by the bilinear part of the generating functional 
S'gg that is proportional to e 2 T 2 , i.e. by the HTL self-energies. 

To construct , or rather its nontrivial part 5S^ 2 \ we sandwich — jj 2 with scalar 
fields and 4i7 /il/ with photon fields. With the leading term of (|2.4|), if rewritten in terms 



of Y, 



one obtains 



5S^ 



-/i 



/3 



(p ip + -m 



HA A + (d AY)^AY 



(A.l) 
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with AY = A - e A and d = = <9 + e V . (A.2) 
Now consider the gauge variations 5 u Aq = d uj and 5^ AY = du . This gives 

<L J I" = J J (lAodou + (d du) - d AY + ud AY^j = , (A.3) 

i.e. gauge invariance of Sjg . The vanishing of (|A.3|) is due to J Q e = in the first term 
(which allows us to rewrite it as 2AYdou) and partial integrations in the two others. 
Therefore, gauge invariance does not require the existence of HTL vertices. Moreover, 
since all HTL's are uniquely determined by f22|, this indeed proves the absence of 
HTL vertices (cf. Sect. 3). 

While the generating functional thus obtained, S e s = S^ s , is gauge invariant, the 
corresponding Lagrangian density is not. A manifestly gauge invariant (mgi) density £' 
can be obtained by adding to £" a suitable total derivative d^J^". It is not hard to construct 
this current as = - J n Y^A^AY. Thus: 

if = e" + d J»= I (d AY -dA )-AY . (A.4) 
Jn d 

( |A.4|) is easily checked to be mgi indeed: S^i' = 0. But in this check J n e = has to be 
exploited, and in fact one can do better. The density £' can be expressed entirely by mgi 
objects, which are the fields E = — VA Q — do A and B = V x A or combinations of them 
as e.g. 

f„ = Y V F VVL = dA^ - d^AY , i.e. (A.5) 
fo = -eE and {f k } = f = -E -e x B . (A.6) 

The mgi action in the QCD case was constructed independently by Frenkel and Taylor 
4~l| and by Braaten and Pisarski f42|| , where the latter just guessed the only possible form 



in accord with known facts and principles. Ref. [flj also gives a proof of its equivalence 
with the non-mgi version. Here we read off the mgi action for scalar electrodynamics from 
eq. (20) of ||2|] by taking its abelian equivalent: 



-W-hf ■ (A7) 



Note that the trace over generators has turned into a factor 1/2. The two versions £' 
and £, both mgi, could differ by a total derivative which is mgi too. However, there is 
no such difference, as we show next. This is in fact already covered by the proof for the 
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non-Abelian case but in the much simpler Abelian case, we can do a straightforward 
calculation which starts with , £' =' and results in '= V . 



In a first step we try to eliminate the potentials from the angular integral £', ( |A.4[ ), in 
favour of the fields E, B . This can be done by using the first two of the following five 
identities: 



d 



XV 



e o e 
d 



-Xd [l + -Z) - 



1 , 
2d' 



where 



X = 



A 



do 
d 



-1 = ixv 1 = 1 

d 2 d ' d 2 d 2 ' 

^ = X(2 + 3Z) + 1(1 + Z) , 



and Z = lv o V - 1 = — V x (V x . . .) 
A A 



(A.8) 
(A.9) 
(A. 10) 

(A.ll) 



( |A.8| ) to ( A.10 ) are replacements allowed under the angular avergage J n if no further 
dependence on e occurs. Using (|A.8|) and, in the last step, Z A = (1/A)V x B one 
obtains 



^ e , -± e o e -? 
E —A — E — - — A 
d d 



E XV A + E Xd 1 + -Z ) A +E —Z A 



1 



2d 



E 



(A.12) 



We get rid of do through the homogeneous Maxwell equation 8qB = —V x E, i.e. 
<9o(l/A)V x B = —ZE. Now ([A.10|) , used in backward direction, and (|A.11[) lead 
to 



2 Jn 



(A.13) 



All three terms are contained also in the desired result £. But two terms are missing. So, 
in the last step, we exploit that 



e x B 



1 

l 2 



e x B + e x B 



d 2 



E = . 



Adding this under J n in ( |A.13|) the calculation ends up with 



1 

2 Jn 



d? 



eE —eE — E +e x B ) — [ E +e x B 



d 2 



(A.14) 



(A.15) 
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see ( [OP with (|A~6|) . 



To verify (|A.14|) , we use (^]9) and ([A. 10 ), note that Z B = — B and eliminate E by 
means of the homogeneous Maxwell equation: 



B^-B-B e^eB-B^r X E 
d z d z d z 



where (|A.9|) explains the last step. ( |A.1|) with 
72|) given in the main text. 



B ^-B +B XB -B X-V x E 

d z d 

^M 1 -?) 5 - 0, (Ai6) 

l" _, f = l and ( |A~7D constitute the result 



Appendix B 



Here the next-to-leading order contributions to the polarization functions U\g and n t 
as given by Eqs. ( p\4p and ( P75| ) are evaluated for all and g. The strategy for doing this 
was already outlined at the beginning of section 6. We start with the imaginary parts of 



the two sums (|6.2| ) and apply formula (|6.4j) . For convenience we write 

S^^A-A/, (i = l,2) with f l = l and f 2 = k 2 - { -^^. (B.l) 

The two integrals in (|6.4j ) allow for the substitutions x uj — x and/or k — > g — A;. 
Doing both transformations the integrand remains unchanged, since the fi are invariants 
under this transformation and the density p is an odd function of x. Thus, the prefactor 
uj may be taken inside the integral as uj = x + (uj — x) — > 2x : 

SmSi = 2ttT f— ) [d i kf i [dx—^—p^{x ) k)p(x-uj ) k). (B.2) 
\2nJ J J x — to 

We insert the spectral densities of ( |6.1|) . After some trivial manipulations with the delta 
functions (with the aim to remove uj_ from the prefactors) we obtain 



T r°° r 1 1 / \ 

SmSj = — / dkk 2 du fi sgn (u; — uj) 6 ylkqu + q 2 — uj 2 + 2ujuj) — (uj — > —uj) 
4tt Jo J-i uj z \ > 

T f°° k 

= - — / dk—^fi sgn (uj - uj) 6 - (w -> -w) , (B.3) 
8nq Jo uj z 
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where, while fi will always equal 1, the function f 2 has turned into 



f2 = k 2 ~ 



4q 2 



q 2 — uj 2 + 2lulu 



(B.4) 



6 stands for the step function O = 9 (2kq — \uj 2 — 2luuj — q 2 \ ). Next we change to the 
integration variable x = \/ fi 2 + k 2 ' = uJ and manipulate the step function: 



= e 



2 2 

q — uj 



^ T 2 



4cj 2 



(B.5) 



where Q 2 = uj 2 {uj 2 — q 2 — 4fi 2 )/(uj 2 — q 2 ) as given in ( p.llj ) in the main text. Obviously 
(and notwithstanding the remaining integration), there is no imaginary part in region II, 
i.e. in q 2 < uj 2 < 4/i 2 + q 2 . For a convenient formulation in the other two regions we 
introduce 



Xl = ^-^n and x 2 = U - + ^n 
2 2u 2 2u 



so that 



6 ( [x — x\ ] [x 2 — x] ) 
O = <j 

( [x — Xi] [x — x 2 ] ) 



in III 

in II 
in I 



(B.6) 



(B.7) 



The function f 2 now reads 



uj 2 - q 2 
q 2 



4/i 2 g 2 + (u z - q 



LUX — X 



2\2 



A(u 2 - q 2 ) 



(B.8) 



and the expression for Sm J2i so f ar obtained is 

T r°°. 1 



2 roo i 

SmSj = / dx — fi sgn (a> — x) — (a; 

87rg 7/i a; 



(B.9) 



Clearly, the integrals ( |B.9| ) can be evaluated analytically, once the range of the parameters 
is specified. Consider region III and assume uj > 0. Then fi < x\ < x 2 < uj. Hence, 
sgn (uj — x) = +1 and the term — (uj — » — uj) vanishes. Thus, in region III: 



5 3m Si 



T 



J 



, uj 2 -q 2 (qVL 4fi 2 q 2 + (uj 2 - q 2 ) 2 ' 
87rg| -^{-2 A(uj 2 -q 2 ) J , 



l = 1 
i = 2 



(B.10) 



with J7" = In | (c<j 2 + qVt)/(uj 2 — qfl) | as given in ( |b.9D in the main text. Of course, in 
region III neither the prefix 5 nor the absolute value in the argument of the logarithm are 
necessary. See however below. 
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In region I and for uj > we have x\ < \x < X2 and sgn {u — x) = — 1. Hence, 
in the first term of (P-9| ), the x-integration runs from X2 to 00. The second term, i.e. 
— (uj — > —00), needs one more position: x\ = — f + <fj^ • Note that < \i < x\ < X2 ■ 
Also, /2 might be split into its even and odd part with respect to uj. In region I and for 
positive uj we obtain 



T 

%itq 



/ dx~f. 

J xi X 



even 
2 



uj 2 — q 2 



■ uj 



dx + dx 

XI J X2 



(B.11; 



Obviously, this expression needs subtraction of the hard leading-order imaginary part, 
which is in fact non-zero in region I. This amounts to subtracting 2 f£°dx from the two 
diverging integrals in (|B.11|) . The result agrees precisely with ( |B . 1 0|) , but the prefix 8 and 
the absolute value in J are no longer superfluous. For i = 1 there is only the first integral 
in ( p. 11|) . For an immediate check, whether the above partial results are correctly used 
in the main text, the reader may take the imaginary parts of ( |3.6|) , ( p.7|) , insert (p.lOj ) 
and compare with ( |6.6| ) and ( |S.5[ ). 



We turn to the real part. It has to be determined from the dispersion relation (|6.3|) . 
But instead of directly treating the corresponding integrals (a torture) there is a pleasant 
way by guessing and using analytical properties. First of all, we put a 5 in front of both 
sides of ( |6.3| ) to indicate subtraction of leading-order terms. The integration over t runs 



from —00 to —y/Afi 2 + q 2 ' , from — q to q and from y/Afi 2 + q 2 ' to 00. This fact may be 
dealt with by including the step function 



e' = 



q 2 -t 2 



4/i 2 + q 2 - t 2 



[B.12) 



Let us start with the case % — 1, i.e. with / = 1 in ( |6.3|) . The imaginary part is an odd 
function of t. To see this explicitly, and by using ( |B.10[ ), (|6.9|), ( |6.11| ), we write ( |6.3| ) as 



8n 2 q 



dt 



t — UJ 



L(t) e' 



with 



L(t) = In 



t 2 + qtW(t) 
t 2 - qtW(t) 



and W(t) 



lA/j 2 + q 2 -t 2 



q 2 -t 2 



(B.13) 



(B.14) 



Principal values are understood whenever a real variable runs over a pole. In order to get 



rid of the absolute value signs in (|B.14|) we rewrite L as 

_ 1 ( 4f*>+(tw+ q r 



(B.15) 
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Next we tackle with the step function in (p.l3| ). Note that dteW(t) vanishes automatically 
in the unwanted region q 2 <t 2 < Aji 2 + q 2 . The same holds true for 8leg(t) with 



9W 



In 



2/i — izW(z) — iq 



(B.16) 



2/i - izW(z) +iq J 

Moreover, in the other regions dteg(t) agrees with ( B.15 ). Thus, in all regions 

L{t)& =9teg{t) = g x . (B.17) 

At this point we recall that 9?e£j and SmEj are originally defined by approaching the real 
axis from the upper half complex plane (UHP). Note that, with e a positive infinitesimal, 
W(t + is) turns into +i|W(t)| when entering region II from larger as well as from lower 
t. Both the functions W(z) and g(z) have cuts on the real axis, but they are analytic in 
the UHP. Since, in addition, g(z) behaves as 1/z at z — > oo, we may state 



dz' 



2-rri Jc z' — z 



g(z') = g(z) 



(B.18) 



where C surrounds the UHP counterclockwise. Through z — > t+ie and g(t+ie) = g±+ig2 
([BT8D tells us the common dispersion relations 



02 (*) 



(B.19) 



Combining the left of these equations with (|B.13|) , ( p,17|) and ( p,16|) , we have 



KeEi = - 



T 
iirq 



g 2 {uJ + ie) 



T 

8^q 



K 



(B.20) 



with 1Z given by (|6.8|) . The procedure is unique. One can now verify 511 g as given in 
II- 

In treating the case i = 2 we follow the same lines of reasoning. <59fte£2 is given by the 



integral (|B.13|) but with L(t) replaced by the function 
M(t) 



t2 - q2 tw(t)- 4 ^ 2+ } t2 - q2)2 m ■ 



2q 



4q 2 



(B.21) 



To restore convergence of (P-13| ) we subtract (and add) on both sides the same expression 
taken at fj, — 0, i.e. we split 



fB.22) 
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with Sfjfte^ = 3?eS2 — KeS^ -0 , and treat the second term of ( |B.22|) at the end. This 
decomposition was used recently in the QCD case (§4 of with the terms called 'one- 
loop soft' and 'one-loop hard', respectively. We observe that 5^M(t) 6' is the real part 
3?e g(u + ie) of the following complex function: 



9{z) 



> g 2 )- [ W(z) - 1 ] - + ( ! 2 " ^ In ^ ~ IZW{Z) ~ * 



2q 



Aq 2 



2/i - izW(z) + iq 



Aq 2 



z + q + ie 
z — q + ie 



2^ 



q 



[B.23) 



where the last term, which obviously is not fixed by gi, has been determined from the 
requirement g(z) ~ 1/z at z — > oo. The relations ( |B.18| ) and ( |B.19| ) hold true again, and 
thus 



8^ Jie E 2 



8n 2 q 



dt- 



t — UJ 



9i(t) 



T 



8nq 



g 2 (uj + ie) 



(B.24) 



Note that uW(u + ie) is equal to the quantity £, ( |6.10| ), used in the main text. With 
g 2 (ui + ie) taken from (|B.23| ) we end up with 



<L3Re£ 2 



T 

8^q 



(B.25) 



2g 1 ' " lr/- 4g 2 

where 0/ = 6(q 2 — uj 2 ) restricts the last term to region I. 

It remains to study the second term of ( B.22|) . Its imaginary part may be read off 
from (IRTOD , 

5^smJ2 A A f 2 



T cu 2 -q 2 



2qu) - (to 2 - q 2 ) In 



8nq Aq 2 

and the appropriate, in the UHP analytic function is 

z + q — ie 



uj + q 




uj — q 





[B.26) 



(Z 2 ~ q 2 



z 



[z 2 - q 2 ) 2 



Aq 2 



In 



z — q — ie 



q z 



(B.27) 



Again the last term (though in conflict with ( B.26|) ) has been added by hand to make 
g(z) convergent at large z. In reality convergence is restored by the Bose function (which 
one could reintroduce in ( p.9|) by 1/x — > n(x)/T). However, such details (if real on 
the real axis) do not influence the result, as we shall take the imaginary part. In fact, 
^smg{uj + ie) = g 2 = —(uj 2 — q 2 ) n 6/ . Consequently, when the second term of (P-22| ) is 
included, the last term in ( P-25Q simply drops out. 



42 



The results (|6.5|) to (|6.7| ) in the main text, which we have also derived by a much more 
laborious direct calculation, are now readily verified. 



Appendix C 



In this Appendix, the ratio Relli/Q 2 is calculated at the light cone u = q, in order 
to verify the statement (|6.20 ) in the main text. Here the decomposition of the Braaten- 



Pisarski scheme into soft and hard loop momenta, where only the former need resumma- 



tion, fails. Instead, the resummed version (|3.7| ) of the longitudinal polarization Hi must 



be used throughout even though the leading contribution will be seen to arise from hard 
internal momenta. Also, the Bose function cannot be expanded. 

The Bose function may be reintroduced by the replacement T/x — > n(x) in the for- 
mulas of the preceding Appendix B. The latest opportunity for doing so is at (|B.9|) : 



1 r 00 

QmEj = - — / dxn{x) [hziw) — h,2\— oj) ] (C.l) 



&nq Jn 

with h 2 (uj) = fi(uj) sgn (u — x) 6 (C.2) 



and fx — 1. /2 and 6 are given by ( |B.8| ) and ( |B.5| ), respectively. The two sums Sj, 
which constitute II i , are defined in (|B . 1| ) . If we are able to grasp the analytical function 
h(z), which has Qmh(u + ie) = h 2 (uj) and is convergent at large z, we may proceed as in 
Appendix B, but this time inside the real x-integration. Indeed, we find 

q 2 — 2xz — 2kq\ k 
q 2 — 2xz + 2kq J ixq 




h{z) = h{z)-\n [ ±. \ fg ■ o ^ ) - — (z 2 - q 2 - 2xz) 6 i>2 , (C.3) 



where k = \/x 2 — /i 2 '. Only in the case i = 2 there is a term to be subtracted for 
convergence. Thus, with ( B.19|) and h\{uj) = 3?e h(u + ie) we obtain the real parts 

1 r°° 

KeS i = / dxn(x)[hx{u)+hx{-u)\ . (C.4) 



One can now form the linear combination (|3.7| ) of the above real parts, with all details 
filled in. The result is conveniently written down with the integration variable k instead 
of x = \/ fi 2 + k 2 ' = uJ : 



&ell e = (u; 2 -q 2 )-—-J dk -nip) 
oTT q Jo OJ 



2 / 'u 2 - q 2 - 2uu - 2kq 

-(2a; — uj) In — — 

V ; \uj 2 - q 2 - 2uu + 2kq 
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-{2U + u) 2 \n 



'oo 2 -q 2 + 2Zooo - 2kq 
K oo 2 — q 2 + 2cooo + 2kq j 



8kq 



(C.5) 



The expression (|C.5|) is the appropriate place for approaching the light cone. After 
dividing both sides by oo 2 — q 2 = Q 2 , (|C.5|) leads to 

„2 



lim 



ir 2 q 2 



J with J 



dk ^ n{uj) 
o oo 



-k + oo In 



oo + k 



00 



(C.6) 



^ Q 2 

To study the high temperature expansion of this integral J, we subtract and add the 
large- A; limit of the integrand and write 



J=Jo + 5J with J 



dk k n(k) 



(2k s 

-1 + 2 In — 



(C.7) 



In the small difference term 5 J one can use n(x) ~ T/x, so that it is readily evaluated 
to be 5J = Tfi7t/2, which leads to the contribution e 2 Tfi/2irq 2 to tftelli/Q 2 . The leading 
term Jq may be given the form 



u 1+e — u 



7T 2 7T 2 2T\ 1 , 

1 In — + 2 lim - / du 

6 3 I fi I e^o e Jo e u -l 



(C.8) 



where the last term brings in derivatives of the Riemann zeta and Gamma functions [43] . 
Equation ( |6.20| ) is thus obtained. Since accurate values of derivatives of the zeta function 
are rarely found in tables, we note that = —0.569 960 993 . . .. 

We now turn to the imaginary part of the resummed polarization Jig at the level of 
( p.5| ). It derives from ( |C1| ), ( |U.2|) directly and is strictly zero in region II. For u 2 < 

+ q 2 , i.e. in the two regions I and II which have the light cone as the common border, 
we find 

^2 



IIe(u) + ie,q) = —Q 9 [q — oj 



2nq 3 



dec tjc 



.to 



00\ f 00 

n\x — n [ x H — 

2) \ 2 



(C.9) 

Here, the effect of resummation is hidden in the lower endpoint Xq of the integration 
interval: 

1 /4/x 2 + q 2 - oo 2 ' 



Q 



q 2 — oo 2 



(CIO) 



Clearly, due to /i ^ 0, Xq comes down from infinity as q becomes larger than u, giving a 
smooth onset of the imaginary part in region I. For oo < q and q — oo e 2 q, 



e 4 T 2 

QmllAoj + ie, q) = — — exp 

07T 



q 



](q - oo) 



(C.11) 
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so that the imaginary part sets in infinitely smoothly in fact. Alarmingly flC.llj) has the 
wrong sign to give rise to weak damping. As it stands, it would give anti-damping as 
soon as the light-cone is traversed into region I. The simple resolution of this startling 
puzzle is that, in region I, QmIIg(uj — ie,q) = —^mll^uj + ie,q). In region II, where 
the next-to-leading order dispersion curve is still undamped, higher-order corrections will 
generate positive damping. The dispersion curve, when it is followed towards and through 
the light-cone, crosses into region I with Qo = uj(q) — ry(<20, so that only the Riemann 
sheet where Qmlli < can be accessed, which corresponds to Landau damping instead 
of anti-damping. Consequently, for (q — to)/ (eT) <C e 2 one has weakly damped plasmon 
excitations with phase velocities just below 1, whereas further down in region I these 
excitations quickly become overdamped. 



Appendix D 



In this Appendix we add some details to the calculations presented in Chap. 7. 



First, as a consistency check we rederive Eqs. (|7.19| , |7T20|) from an equivalent Schwinger- 
Dyson equation, where the diagrams of Fig. 7 are flipped so that the bare vertices are 
on the right side. Since II is an asymmetric combination of the real-time components, 
this leads to a different starting point for n^J R . In place of Eq. (|7.14|) we obtain 



njS R ( Q ) = n^(Q)+nj$"(Q) 

= \ J y0{ [rf> + ] A R (Q + K)A R (K) 
+ [r^+r^l A A (Q + K)A A (K) 



+ r£« [A R (Q + K)A P (K) + A P (Q + K)A A (K)]}(2K + Q) l/ , (D.l) 
where we have additionally introduced 

a, b 



0,6 



The Ward identities that allowed us to solve T R ^\ T R (°\ and in terms of the 
corresponding self-energies can be derived by subjecting the type-1 and type-2 field to 
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the same gauge transformations [3^]. This would not determine the vertices appearing 
above. Assuming, however, that the theory is separately invariant under type-1 and 
type-2 gauge transformations relates all components of the vertices to the corresponding 
self-energy components. In analogy to (|7.17 ) we are then led to 



T«W(K,K + Q) = -ie(2K + Q)^ 
T^\K,K + Q) = - ie (2K + Q), 
r?°\K,K + Q) = ie(2K + Q) 



1 



2KQ + Q 2 



[S A (K)-~ R (K + Q)\ ,(D.3) 



2KQ + Q 2 
1 

^2KQ + Q 2 " 



S P (K + Q) 
p(K) • 



(D.4) 
(D.5) 



Inserting these into ( p.l|) and keeping only terms involving the Bose distribution function 
exactly reproduces ( |7.20|) . 

In the following we shall concentrate on 77^ in the form ( |7.19| ), since there the con- 
tributions from the poles of the Bose distribution functions that are proportional to T 2 
or to T cancel upon integration of the 3-momenta. Closing the contour around the upper 
half plane, only the poles of Ga{K) contribute, which are located at K = i( ± uo k with 
u k = VF + ^-C 2 '- 

We consider first n QO (Q Q ,0), which when nonvanishing signals the loss of gauge in- 
variance. Always omitting temperature-independent contributions, Eq. ( [7.1 9|) leads to 



#oo(Qo,0) 



dkk 2 



+ 



Q + 2< J 7T 2 
Qo^k 

K 

Qo^k 



n(u k + <) - n(uj k - i() 
(u k + Qo) [n(uj k +i( + Qo) - n(u k + i() ] 
(uj k Qo) [n(uj k Qo) - n(u k - i() ] | . (D.6) 



Taylor-expanding the distribution functions and keeping only terms up to ( 3 , we have 

e 2 i( rdkk 2 



#oo(Qo,0) 



Qo + 2i( 



uj k + i( + Qq 
u k 

u k -K-Qq 



7T 



E 

oo 

E 



m=l . 



2n'(uk) - ( 2 n"'(uj k ) 

n (m) + i£ n (m+l) 

n {m) - i(n {m+1) 



t 2 


Qo- 1 


2 


ml 


C 2 


(-Q ) 


2 


ml 



\m—l 



■(D.7) 
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The integrals associated with higher powers of ( are increasingly infrared-dominated, so 
that the terms neglected are down by powers of (/fi ~ eln(l/e). Collecting the various 
powers in Q , we find for Q «T 



#oo(Qo,0) 



U( 3 e 2 T [dkk 2 



Qo + 2i(J 7T 2 

+ 



nQl m - 2 , (n + l)Q 2m 



(fc 2 + /i 



2W+1 



2W+2 



+ 0(C//x) 



(D.8) 



This sum is telescoping — all terms but the first cancel each other, and ( |7.22|) is readily 
obtained. 

Considering next 77^(0, g), we have 



77 00 (0,g) 



7T- 



3?e / dkk' 



1 



k 2 + /i 2 - 2C 2 + 2iCw fc . 2kq - 4<cu fe + q 2 + 4( 2 
H : m 



n u (0,q) = - [ dkk 2 

TV 2 J 



kq ' ~ -2kq - 4i(u; k + q 2 + 4( 2 



Ak 2 + 8i(uj k - q 2 - 8C 2 , 2kq - AiCu k + q 2 + 4C 2 
In 



4fcg 



-2A;g - 4<a; fc + q 2 + 4( 2 



In the limit q — > this simplifies to 



2 

^oo(0,g -> 0) = -^rKe / 



,m(tjfc + <) 

c 



7T^ 



77^(0,9-^0) = 
which yields Eqs. ( fT23Q . 



e 

7T 2 



- — / Cifc k' 



'C 2 



3 k 2 + fi 2 



n"'(Jk 2 + fi 2 ) + ... 



(D.9) 



(D.10) 



(D.ll) 
(D.12) 
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Z (q) 
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